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Abstract 



The D0 collaboration has performed a study of spin correlation in tt produc- 
tion for the process tt — > bW + bW~ , where the W bosons decay to ev or fiu. 
A sample of six events was collected during an exposure of the D0 detector 
to an integrated luminosity of approximately 125 pb" 1 of yfs = 1.8 TeV pp 
collisions. The standard model (SM) predicts that the short lifetime of the 
top quark ensures the transmission of any spin information at production to 
the ti decay products. The degree of spin correlation is characterized by a 
correlation coefficient k. We find that k > —0.25 at the 68% confidence level, 
in agreement with the SM prediction of n = 0.88. 
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Pair production of top quarks has been observed |J in pp collisions at y/s = 1.8 TeV 
by both the CDF and D0 collaborations, and the mass and production cross section have 
been measured in various channels PH- The observed properties agree well with predictions 
from the standard model (SM). 

For a top quark mass of m t = 175 GeV, the width of the top quark in the SM is 
r t = 1.4 GeV while the typical hadronization scale is Aq C d ~ 0.22 GeV ||. The time 
scale needed for depolarization of the top-quark spin is of the order m t /AQ CD ^> 1/T t 
implying that polarization information should be transmitted fully to the decay products of 
the top quark. That is, the expected lifetime of the top quark is sufficiently short to prevent 
long distance effects (e.g. fragmentation) from affecting the ti spin configurations, which 
are determined by the short distance dynamics of QCD at production |7HlT|. 



The observation of spin correlation in the decay products of ti systems is interesting 
for several reasons. First, it provides a probe of a quark that is almost free of confinement 
effects. Second, since the lifetime of the top quark is proportional to the Kobayashi-Maskawa 
matrix element |Vtb| 2 , an observation of spin correlation would yield information about the 
lower limit on |V^|, without assuming that there are three generations of quark families |L2 



Finally, many scenarios beyond the standard model [p^-[lq| predict different production and 
decay dynamics of the top quark, any of which could affect the observed spin correlation. 

In the decay of a polarized top quark, charged leptons or quarks of weak isospin — \ are 
most sensitive to the initial polarization. Their angular distribution in the rest frame of 
the top quark is given by 1 + cos 9, where 9 is the angle between the polarization direction 
and the line of flight of the charged lepton or down-type quark. Due to the experimental 
difficulties of identifying jets initiated by a down-type quark, we only consider top-quark 
events in dilepton channels, i.e., where both W bosons in an event decay leptonically [ev or 
[iv). The advantages associated with using these channels are that: (1) objects sensitive to 
the polarization of the top quark are clearly identified, (2) background is small compared 
to the lepton+jets channels, and (3) there are fewer ambiguities associated with assigning 
objects observed in the detector to their originating quarks. The disadvantages are that the 
number of events in the dilepton channels is small, and that it is necessary to reconstruct 
two neutrinos in an event whose combined transverse momenta gives rise to the observed 
transverse momentum imbalance in the event. 

At yfs = 1.8 TeV, 90% of the top quark pairs arise from qq annihilation, and, for 
unpolarized incident particles, the produced t and i are also expected to be unpolarized. 
However, their spins are expected to have strong correlation |T2|JT^] event by event and 
point along the same axis in the ti rest frame In an optimized spin quantization basis 



called the "off-diagonal" basis, contributions from opposite spin projections for top quark 



pairs arising from qq annihilations are suppressed at the tree-level |T8[ and only like spin 
configurations survive. This spin quantization basis can be specified using the velocity (3* 
and the scattering angle 9* of the top quark with respect to the center-of-mass frame of the 
incoming partons. The direction of the off-diagonal basis forms an angle ip with respect to 
the pp beam axis that is given by P 



, p* 2 sin 9* cos 9* 
1 — p* z sin 9* 

This particular choice of basis is optimal in the sense that top quarks produced from qq will 
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have their spins fully aligned along this basis. In the limit of top quark production at rest 
{(3* = 0), the t quark and the t quark will have their spins pointing in the same direction 
along if) = 0. 

Defining 6 + as the angle between one of the charged leptons and the axis of quantization 
in the rest frame of its parent top quark, and similarly defining for the other charged 
lepton, the spin correlation can be expressed as |18|,^9|: 



1 d 2 a 1 + k cos 9 + ■ cos 6L 

ad(cosfl+)d(cos#_) = 4 ' ^ 

where the correlation coefficient k describes the degree of correlation present prior to imposi- 
tion of selection criteria or effects of detector resolutions. For tt production at the Tevatron, 



the SM predicts k = 0.88 [21|. In the off-diagonal basis, the correlation coefficient for qq — > tt 
is k = 1. When gg — > tt is included at yfs = 1.8 TeV, the correlation is reduced to k — 0.88. 
The distribution is symmetric with respect to the exchange of 9 + and 8-, and it is therefore 
not necessary to identify the electric charge of the leptons. The physical meaning of k in 
any spin quantization basis corresponds to the fractional difference between the number in 
which the top-quark spins are aligned and the number of events in which they have opposite 
directions. 

The events used in this analysis are identical to those used to extract the mass of the 
top quark in our dilepton sample . They were recorded using the D0 detector p2| , which 
consists of a non- magnetic tracking system including a transition radiation detector (TRD), 
a liquid- argon/uranium calorimeter segmented in depth into several electromagnetic (EM) 
and hadronic layers, and an outer toroidal muon spectrometer. The final sample consists of 
three e\i events, two ee events, and one event, with expected backgrounds of 0.21 ± 0.16, 
0.47 ± 0.09, and 0.73 ± 0.25 events, respectively 0. 

To study the distribution in (cos#+, cos6L), we must deduce the momenta of the two 
neutrinos. The weighting scheme we use is the previously-developed neutrino weighting 
method ||. In this method, each neutrino rapidity is selected from a range of values following 
a distribution consistent with the decay kinematics in tt events. We assume the tt dilepton 
decay hypothesis, and the constraints that m{l\V}) = m(7 2 ^ 2 ) = mw an d m(7iz/i&i) = 
m(7 2 ^ 2 & 2 ) = m t . The problem can be solved by providing a specific input mass m t that we 
assume to be m t = 175 GeV. We then solve for the neutrino momentum vectors, obtaining 
up to four solutions, and assign a weight to each solution to characterize how likely it is to 
represent tt production. A weight is assigned to each solution based on the extent to which 
the sum of transverse momentum components Y^Pk(vv) (k = x, y) of the two neutrinos in the 
solution agrees with the measured missing transverse momentum component ^ (k = x, y) 
in the event. A Gaussian distribution with a width of 4 GeV is assumed for each component 
of the $k [H • The weight is calculated as: 



w 



n ex p 



=x,y 



(gfc -pk{vv)f 
2a 2 



(3) 



The physical objects in the events are smeared to take into consideration the finite resolution 
of the detector, and we consider both possible pairings of the two charged leptons with the 
two jets assigned to b quarks. The presence of a third jet is also taken into consideration [H. 
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TABLE I. Asymmetry values for the 6 dilepton events at D0 



Event Number Event type A 

10822 ee 034 

12814 efi -0.16 

15530 fj.fi 0.50 

26920 efi 0.85 

30317 ee 0.52 

417 efi -0.19 
< A > 0.31 ±0.22 



For each solution, we can then boost the decay products into the rest frame of the 
original top quarks and calculate the relevant decay angles (cos 9 + , cos 9-). The event fitter 
returns many such solutions for an event, and the goal is to deduce the original value of 
(cos#+, cos6L) from the reconstructed distributions. 

The differential cross section depends on the product £ = cos#+ ■ cos6L. We define an 
asymmetry A for all solutions in an event as |7],|l(],|l^] : 

A= i(f 1 ^.r^A. (4) 



For perfect resolution and acceptance, A is expected to be «/4. 

Since the event fitter returns solutions with assigned weights and there is no "unique" 
solution, we sum the weights for all the solutions to populate the distribution £, which is 
shown in Fig. [l| for the 6 events. The values of A are listed in Table |I[ 

Monte Carlo event generators such as herwig [23] and pythia ]24J, in their current 
implementation, do not take proper account of spin correlation in ti production, and the 
two top quarks in an event are made to decay independently of each other, i.e. k = 
is assumed. To include the effects of spin correlation, ti events from the pythia event 
generator are sampled at the generator level with the weight (1 + /«£), where £ is calculated 
from information at the generator level. We have checked this method against a Monte 
Carlo containing a fully correlated matrix element (where k — 1 for ti events initiated from 
qq annihilation) and found the two methods are equivalent ||21|| . 

To estimate the sensitivity of our method, we created 1500 ensembles of 6 events con- 
sisting of appropriate fractions of ti signal and background. From Monte Carlo studies, 
we expect A = 0.207 ± 0.006 for full spin correlation (k = 1) when all detector and back- 
ground effects are included, while A = 0.25 for perfectly reconstructed events without any 
background. The statistical uncertainty on our measurements is estimated to be 0.20 from 
these ensemble studies. Similar tests were performed for ensembles of 6 events without spin 
correlation (k = 0), and we find an expected A = 0.115 ± 0.005, while ideally ^4 = 0. The 
main cause for loss of sensitivity is the incorrect pairing of the lepton with the jet. This 
produces a strong bias in A [ET[. From the Monte Carlo samples generated with values of k 
between —1 and 1, we find a linear relationship between A and k: A = 0.112 + 0.088k. 

We obtain A = 0.31 ± 0.22 from our data, which translates into k = 2.3 ± 2.5, assuming 
that a linear relationship between A and k also holds beyond — 1 < k < 1, though the values 
\k\ > 1 are not physical. 
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FIG. 1. Distribution of £ for the 6 dilepton events. 



Systematic uncertainties are negligible compared to the statistical uncertainty in our 
result. Varying the top quark mass by 5 GeV results in a shift in A of 0.01. There has 
been no theoretical calculation of effects of gluon radiation on the spin correlation of the 
top quarks. However, these effects were studied for spin-uncorrelated events (i.e. n = 0) 
by including gluon radiation in the pythia event generator. This results in a shift in A 
of 0.0065 ± 0.0063, where the error is due to finite Monte Carlo statistics. The asymmetry 
distribution expected from background is similar to that for spin-uncorrelated ti events, and 
its impact is small. 

To maximize the physical information present in the data, the full two-dimensional phase 
space of (cos#_|_, cos#_) is used in a two-dimensional binned likelihood analysis. The phase 
space is split into a 3 x 3 grid, each side of which spans 1/3 of the range of cos^ and 
cos6L. The nine bins are populated for data with weights (wi, ...,wg) from the event fitter, 
with the distribution of weights for each event normalized to unity. Similar distributions 
are made for the generated Monte Carlo events using different values of k for ti signal and 
an appropriate admixture of background. Comparisons of data with Monte Carlo are used 
to extract n. 

Because an event populates each bin with fractional probability, a simple likelihood 
assuming a Poisson distribution may not be appropriate. Moreover, since the weights for 
each event satisfy the normalization condition = 1, only eight out of the nine weights 

are independent, and there are correlations among the weights in any given event. 



To find eight independent variables, the covariance matrix Cij 



COV [Wi, w j/ 
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FIG. 2. Plots of probability density for tt events in the dilepton channels in (cos 8+, cos 
phase space. Top left: Monte Carlo events with k = —1; top right: Monte Carlo events with 
K = +1; bottom left: our data; and bottom right: the likelihood as a function of K showing the 
68% confidence limit of k > —0.25. The box area is proportional to the summed weights in the 
bin. 

1, . . . , 8) is calculated from the Monte Carlo events for a given spin correlation k and 
background, and diagonalized using a matrix A, such that A~ X CA has only diagonal el- 
ements. The new independent variables (i.e. diagonalized weights) are found by applying 
this transformation matrix to the weights, V = A~ X W ', where W = (wi, . . . , u>s) T an d 
V = (v i, . . . , v 8 ) T . The distributions (i = 1, . . . , 8) of the new variables Vi are used to 
define the likelihood 

N 8 

£(*)=nn ( 5 ) 

i j=l 

where % are the new variables for ith event and is the number of events. By explicitly 
constructing the likelihood, we do not have to make any assumptions about the underlying 
distributions of the weights. 

The result is shown in Fig. [|. The probability densities for the Monte Carlo generator 
at k = — 1 and k — 1 are shown for comparison. From the dependence of the likelihood on 
k, we can set a 68% confidence interval at k > —0.25, based on the line fit, in agreement 
with the SM prediction of k = 0.88. 

In conclusion, we have presented a search for spin correlation effects in the production 
of tt pairs in pp collisions at ^/s = 1.8 TeV, where the dominant production mechanism is 
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expected to be the annihilation of incident qq states. This analysis makes use of the fact 
that there exists an optimal spin quantization basis for the produced top quarks, and that 
the charged leptons from top-quark decays are most sensitive to the polarization of the top 
quark. From this analysis, we conclude that k > —0.25 at the 68% confidence level, which 
is compatible with correlation of spins expected on the basis of the standard model. 
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